Computing a Matrix Function for Exponential Integrators
نویسنده
چکیده
An efficient numerical method is developed for evaluating φ(A), where A is a symmetric matrix and φ is the function defined by φ(x) = (ex − 1)/x = 1+ x/2 + x2/6+ .... This matrix function is useful in the so-called exponential integrators for differential equations. In particular, it is related to the exact solution of the ODE system dy/dt = Ay + b, where A and b are t-independent. Our method avoids the eigenvalue decomposition of the matrix A and it requires about 10n3/3 operations for a general symmetric n×n matrix. When the matrix is tridiagonal, the required number of operations is only O(n2) and it can be further reduced to O(n) if only a column of the matrix function is needed. These efficient schemes for tridiagonal matrices are particularly useful when the Lanczos method is used to calculate the product of this matrix function (for a large symmetric matrix) with a given vector.
منابع مشابه
A Krylov Subspace Algorithm for Evaluating the Φ-functions Appearing in Exponential Integrators
We develop an algorithm for computing the solution of a large system of linear ordinary differential equations (ODEs) with polynomial inhomogeneity. This is equivalent to computing the action of a certain matrix function on the vector representing the initial condition. The matrix function is a linear combination of the matrix exponential and other functions related to the exponential (the so-c...
متن کاملNorges Teknisk-naturvitenskapelige Universitet Expint — a Matlab 1 Package for Exponential Integrators
Recently, a great deal of attention has been focused on the construction of exponential integrators for semi-linear problems. In this paper we describe a matlab package which aims to facilitate the quick deployment and testing of exponential integrators, of Runge–Kutta, multistep and general linear type. A large number of integrators are included in this package along with several well known ex...
متن کاملA matlab 1 package for exponential integrators
Recently, a great deal of attention has been focused on the construction of exponential integrators for semi-linear problems. In this paper we describe a matlab package which aims to facilitate the quick deployment and testing of exponential integrators, of Runge–Kutta, multistep and general linear type. A large number of integrators are included in this package along with several well known ex...
متن کاملThe Leja Method Revisited: Backward Error Analysis for the Matrix Exponential
The Leja method is a polynomial interpolation procedure that can be used to compute matrix functions. In particular, computing the action of the matrix exponential on a given vector is a typical application. This quantity is required, e.g., in exponential integrators. The Leja method essentially depends on three parameters: the scaling parameter, the location of the interpolation points, and th...
متن کامل